home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlasd5.z / dlasd5
Encoding:
Text File  |  2002-10-03  |  3.7 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSSDDDD5555((((3333SSSS))))                                                          DDDDLLLLAAAASSSSDDDD5555((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLASD5 - subroutine computes the square root of the I-th eigenvalue of a
  10.      positive symmetric rank-one modification of a 2-by-2 diagonal matrix
  11.      diag( D ) * diag( D ) + RHO * Z * transpose(Z)
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DLASD5( I, D, Z, DELTA, RHO, DSIGMA, WORK )
  15.  
  16.          INTEGER        I
  17.  
  18.          DOUBLE         PRECISION DSIGMA, RHO
  19.  
  20.          DOUBLE         PRECISION D( 2 ), DELTA( 2 ), WORK( 2 ), Z( 2 )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      This subroutine computes the square root of the I-th eigenvalue of a
  37.      positive symmetric rank-one modification of a 2-by-2 diagonal matrix
  38.      diag( D ) * diag( D ) + RHO * Z * transpose(Z) . The diagonal entries in
  39.      the array D are assumed to satisfy
  40.  
  41.                 0 <= D(i) < D(j)  for  i < j .
  42.  
  43.      We also assume RHO > 0 and that the Euclidean norm of the vector Z is
  44.      one.
  45.  
  46.  
  47. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  48.      I      (input) INTEGER
  49.             The index of the eigenvalue to be computed.  I = 1 or I = 2.
  50.  
  51.      D      (input) DOUBLE PRECISION array, dimension ( 2 )
  52.             The original eigenvalues.  We assume 0 <= D(1) < D(2).
  53.  
  54.      Z      (input) DOUBLE PRECISION array, dimension ( 2 )
  55.             The components of the updating vector.
  56.  
  57.      DELTA  (output) DOUBLE PRECISION array, dimension ( 2 )
  58.             Contains (D(j) - lambda_I) in its  j-th component.  The vector
  59.             DELTA contains the information necessary to construct the
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAASSSSDDDD5555((((3333SSSS))))                                                          DDDDLLLLAAAASSSSDDDD5555((((3333SSSS))))
  71.  
  72.  
  73.  
  74.             eigenvectors.
  75.  
  76.      RHO    (input) DOUBLE PRECISION
  77.             The scalar in the symmetric updating formula.
  78.  
  79.             DSIGMA (output) DOUBLE PRECISION The computed lambda_I, the I-th
  80.             updated eigenvalue.
  81.  
  82.      WORK   (workspace) DOUBLE PRECISION array, dimension ( 2 )
  83.             WORK contains (D(j) + sigma_I) in its  j-th component.
  84.  
  85. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  86.      Based on contributions by
  87.         Ren-Cang Li, Computer Science Division, University of California
  88.         at Berkeley, USA
  89.  
  90.  
  91. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  92.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  93.  
  94.      This man page is available only online.
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.